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Abstract 

Clusters of galaxies are massive structures containing approximately 10^^ solar 
masses. The mass contained in clusters provided the first hints of the existence of 
dark matter. The hot intra-clustcr medium (ICM), which is a gas with a temperature 
of 10® Kelvin, is easily detected by X-Ray telescopes. Of particular interest is the close 
Perseus galaxy cluster, one of the brightest X-Ray emitting galaxy clusters in the sky. 
Due to the large angular extent of the cluster, background noise subtraction is particu- 
larly challenging, and has yet to be completed using the full set of Perseus observations 
from the XMM-Newton X-Ray Telescope. Detailed temperature and abundance radial 
profile maps have revealed a significant lack of homogeneity within the cluster. Previ- 
ous surveys of Perseus with the Suzaku telescope, which has a worse angular resolution 
and less light collecting area than XMM-Newton, revealed over-densities of X-Ray 
emission. These results provide evidence that the baryon fraction exceeds the univer- 
sal average, which we had initially hoped to study. We have yet to confirm or deny 
the existence of clumping in these regions, which could explain such over-abundance 
of X-Ray emission. 

This project offers a framework of efficient, automated processing techniques to 
"clean" images of noise from the mechanics of the telescope, background radiation from 
local sources such as the solar wind, and more distant sources such as background AGN. 
The galaxy cluster studied in this project contains high levels of contamination due to 
its line-of-sight position close to the dust- and star-filled arms of the Milky Way galaxy. 
Rigorous spectral model fitting of the cluster employ multiple parameters dedicated 
to accounting for these contaminations. The framework created from this analysis 
technique will provide the opportunity to expand this analysis to any nearby galaxy 
cluster, such as the Virgo, Coma, and Ophiuchus Clusters. This research should provide 
significant insight into how matter, both baryonic and dark matter, is distributed 
throughout diffuse cluster systems, as well as give clues to the origin of the ICM. 
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1 Introduction 



1.1 Galaxy Clusters 

The study of X-Ray images of galaxy clusters is a relatively new field of astrophysical 
study that has benefitted greatly in the past decade from the launch of space based X-Ray 
observatories. Galaxy clusters are large, gravitationally bound systems consisting of many 
galaxies. The study of the physical characteristics of galaxy clusters offers key insight into 
cosmological parameters, such as the distribution of dark matter in the universe, the effects 
of the cosmological constant, and the distribution of baryonic matter in these enormous 
structures. Clusters often contain a central active galactic nuclei (AGN), which is theoret- 
ically powered by a matter-accreting supermassive black hole. Clusters themselves can be 
active as well; recently having experienced a collision with another cluster or adding a new 
component galaxy could be a potential cause for the cluster ICM to not be in hydrostatic 
equilibrium. Clusters that are in a relaxed hydrostatic equilibrium ([T]) state are typically 
modeled as spherically symmetric systems. 



Clusters offer numerous potential avenues of study. Clusters are used as model laboratories 
for investigating large scale creation and evolution of elements, as the large amounts of 
interacting gas are theorized to originate from stellar fusion processes. The extremely thin 
quantity of hot gas (< 10~^ particles/cm~^) in the ICM has emissions well within the X-Ray 
wavelength range. The composition and exact origin of this gas is an active area of research, 
which this study aims to explore in depth. 

1.2 Cluster Composition 

The main elemental composition of the ICM is important for a variety of reasons. The 
exact abundance of specific elements, such as iron, provides researchers with the ability to 
model the age and necessary mass of progenitor stars. Iron, which is produced as the upper 
fusion limit in high-mass stellar cores, can also be produced by the violent explosion of white 
dwarf stars in Type la supernovae. Mergers of the cluster with new component galaxies 
infuse the colliding ICM with high levels of supernova ejecta. By measuring the distribution 
of iron abundances, insight can be made into the evolutionary past of the star forming 
regions of the component galaxies of the cluster [2]. Similarly, alpha element abundances, 
which consist of O, Mg, Si, S, Ca and Ti, are of great importance to the study of the origin 
of the first stars and supernova. Following the more recent evolutionary track and formation 
of elements provides a more complete picture into how the current elemental composition 
of the universe came to be [3]. Recent studies have indicated that the majority of the iron 
abundance observed is directly formed in Type II supernovae, models of which are used to 
provide constraints on the rate of present and past star formation. This places limits on the 
timescale during which the iron and alpha element enrichment could have occurred |8j. 
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1.3 Mass Determination 



A frequently used technique in mass determinations is to assume the ICM of the cluster in 
question exhibits hydrostatic equilibrium as well as spherical symmetry. Using the model of 
hydrostatic equilibrium, one can derive a homogenous fit of gravitating mass across the ob- 
served portions of the cluster [TT]. Such mass determinations have previously been emission 
model-dependent, relying upon multiple a priori assumptions concerning the state of the 
cluster beyond the aforementioned hypothesis of hydrostatic equilibrium and spherical sym- 
metry. Recent advancements have been made to eliminate some of the a priori assumptions 
made in the case of mass determinations. In particular, the Navarro, Frenk, and White mass 
density model revealed that the density of dark matter halos around galaxy clusters, which 
form during the relaxation stage following mergers and collisions, do not necessarily depend 
on cosmological constants (A and flm) or the power spectrum fit applied to the dark matter 
halo [9]. In recent years, even more significant advancements have been made, culminating 
in the development of cluster mass mixing models for X-Ray observations. In particular, the 
CLMASS and NFWMASS models provide a robust method of producing density and gravitational 
potentials for the spherical shells of clusters in hydrostatic equilibrium. In particular, when 
high quality observational data is available, these models provide significant improvements 
over the classical approach to mass determination [TO] . 

1.4 Motivation for Study 

This study concerns the analysis of archival data from the Perseus galaxy cluster. In a 
recent study by Simionescu et al., Suzaku data were used in conjunction with the CLMASS 
and NFWMASS models. The results from these fits reveal that the baryon gas fraction of the 
outer regions of the cluster (exceeding 2/3 of the virial radius) exceeds the cosmic average 
baryon gas fraction [13j. Such a result confiicts with accepted theory and provided additional 
motivation for this study of the Perseus Cluster. 

The archival data obtained for this study is from the XMM-Newton telescope. This 
telescope, in comparison to Suzaku, offers significantly better angular resolution due to 
XMM-Newton's better optics. Suzaku had two main benefits, a shorter focal length and a 
lower radiation environment. These advantages give Suzaku a lower background emission 
profile. Suzaku also has a smaller light collecting area compared to XMM. With a higher 
angular resolution, evidence of clumping will become significantly more apparent when ap- 
plying mass determination and density models to derived spectra. The larger effective light 
collecting area also provides the benefit of higher photon counts over a shorter exposure 
length, providing better overall statistics. This should offer the capability to detect any 
clumping of matter in regions noted to have higher-than-anticipated baryon gas fraction. 
If such clumping effects are observed, this could easily explain the conflict present in the 
Simionescu et al. study. Through this study, it is the goal of the investigator to supply am- 
ple reason for the continued study of diffuse and extended galaxy clusters using the analysis 
procedures to follow. 
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1.5 Expected Issues 



The study of Perseus using the XMM-Newton observatory poses several unique problems. 
The field of view of the instrumentation on board the telescope is roughly 30'. The cluster 
itself has an extent of nearly 2°, much larger than the FOV of the telescope. Because of the 
extent of the cluster, there is no region that can easily be selected to adequately model the 
background emission of the cluster. This issue is dealt with through the use of the ESAS 
data reduction package. This recently released analysis package provides several new tools 
to deal with background effects from extended diffuse X-Ray sources. This involves multiple 
tools designed particularly for the purpose of cleaning processed event files of soft proton, 
quiescent particle background, solar wind background, and crosstalk contamination (due to 
bright X-Ray sources located off-axis from the central observation heading). Furthermore, 
the novel use of ROSAT All-Sky Survey data as a spectral model to diffuse X-Ray background 
provides a powerful method to eliminate background emissions. Due to the very recent release 
of both the CLMASS and ESAS packages, this method of analysis has yet to be attempted on 
the Perseus galaxy cluster. 



2 Methodology 

2.1 Initial Data Reduction 

The data processing techiques used to create spectra and images of the cluster are based 
off of the XMM-Newton Extended Source Analysis Software (ESAS) [H]. Significant im- 
provements were made to the workfiow of this software distribution. In particular, a fully 
automated framework was designed, based on work done in [16], to streamline data pro- 
cessing. Mechanisms to catch and eliminate errors were also introduced. This is discussed 



in more detail in ^2.5 



The ESAS software package first examines and produces a set of unprocessed and un- 
cleaned images directly from the detector with the odf ingest routine. This is followed by 
the production of calibration files by the cifbuild routine, which allows reference calibra- 
tion detector files to be used in the production of cluster and background spectra as well as 
image files. 

A basic preparation task, emchain, first generates the event lists. Then, mosf ilter (an- 
other utility), examines the unprocessed files for signs of obvious soft-proton contamination. 
It also provides a master clean file that is used to derive spectra. The soft-proton con- 
tamination characteristics can be examined directly through the plotting of the observation 
histogram and light curves (see Fig. [T|. Observations that contain low levels of soft-proton 
events tend to have a light curve distribution that is generally uniform. Observations that do 
experience significant soft-proton contamination tend to observe low levels of high intensity 
events. Such contaminated observations in this study were noted, as further spectra fitting 
would indicate whether the not reliable results could be produced, or if the observation under 
study had to be discarded. 
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Figure 1: Light curves of the MOSl detector for the first exposure of ObsID 0085110101. 



2.2 Masking and CCD Examination 

The first step to producing spectra involved the creation of point source mapping. As 
the primary processing techniques emphasized in ESAS are for diffuse spectra, point sources 
can cause a significant level of contamination in the resulting spectra. The ESAS cheese 
task was used to exclude such sources. The cheese utility runs source-detection algorithms 
to locate point sources, which are then added to a point source list. From this list, image 
masks (see Fig. [2]) are created to exclude these point sources from final images. 

Point sources were defined by the parameter selection of cheese. The threshold intensity 
value of an 10~^^ ergs/s cm^ was used for the detection of potential point sources. Detected 
sources were then excised by removing the region radially surrounding the source to half 
the surface brightness of the background. The source detection algorithm was limited so 
that sources within 40" of each other were not excluded, so as to prevent brighter, extended 
sources from excision. However, this required the manual examination of the image for 
bright, closely-packed foreground sources, which the cheese task would skip. Any such 
sources were then excluded. 

For the purpose of this study, images were generated by cheese and all subsequent tasks. 
Masks generated by cheese were made for the range of energies between 400 eV through 10 
keV. 

Detector images of the observation were then examined for CCDs operating in anoma- 
lous states. CCDs behaving in this format tend to over-estimate low energy (< 1 keV) 
background. As the Perseus cluster has significant low energy emission (~ 0.3 keV), ex- 
clusion of anomalous CCDs states was vital. CCDs were determined to be in anomalous 
states through the output of the corner hardness ratios in the mos-f ilter task as well as 



6 



examining the images visually. A shell script was used to extract the notably bad CCDs 
from the output of mos-f ilter and automatically set future tasks to exclude them from 
processing. 




Figure 2: Cheese mask image of ObsID 0085110101. This image indicates the overall profile 
of the MOS detectors. 



2.3 Annulus Generation 

Production of the annular shaped sections was necessary for CLMASS and NFWMASS mixing 
model analysis. This required the creation of several Python-based programs to create 
spectra across identical sky coordinates for both MOS detectors. Instead of creating masking 
regions for the included mos-spectra and mos_back tasks, partial annuli of events were 
computed directly from event files. This involved use of the kapteyn Python module and 
a novel rotation algorithm described below. An example of the resulting event files can be 
seen in Fig. |3j When generating an annular section, the events contained by the annulus 
section were rewritten into a new event file containing the identical FITS formatted header 
of the original event file. 

This program can be run in two modes, generating events of full cluster annuli or annuli 
sections propagating in a specific radial direction from the cluster center. This system 
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Figure 3: Left: A partial annulus from a non-center observation. Right: A full annulus from 
the center cluster observation. 



permits a wide variety of potential analysis and lends itself particularly well to the use 
of spectra in mass mixing models (such as CLMASS and NFWMASS), which determine the 
characteristics of the spherical shells defined by the cylinder of three-dimensional emission 
described by the annuli (Fig. |4]). 

Event files that were generated were sorted according to the radial spacing between 
the great arcs describing the annulus as well as the angle in radians described by the flat 
edges of the annulus. This is particularly useful when generating spectra of the full annulus 
encompassing the cluster, as it determines the exact angular extent of each annulus. In the 
case of directional analysis, this is particularly beneficial for any situation where a direction 
chosen might lie along the edge of the observations field of view; the angle is defined by the 
maximum and minimum right ascension of all events within the reformatted event file. This 
angle generation is vital for the use of many mixing models, as it is necessary to describe 
the particular region over which to apply the model. 

2.3.1 Annulus Section Selection Algorithms 

Two annulus selection algorithms were used in this research: 

• partialAnnulus: Used to compute an annulus section S given an inner radius ri, an 
outer radius r2, and two bounding angles 6i and 62- 

• fullAnnulus: Used to compute an annulus A given by radii ri and r2, along with an 
associated bounding angles in the case that the annulus was not complete (i.e. when 
the annulus was generated from a peripheral, non-center observation). 



Figure 4: Diagram of spherical shells of cluster from [TO]. 



partialAnnulus{evtfile, Pq, Sq, ri, r2, Oi, 62): 

rotevts ^ map rot{po, So,-,-) on evtfile 

return {p', 6' G rotevts '■ 0i < p' < O2 Sz ri < 6' < r2} 

Table 1: A pseudo-code implementation of the partialAnnuli routine in mathematical short- 
hand. 

fullAnnulus{evtfile, po, 60, ri, r2): 

rotevts map rot{pQ,6o, -, ■) on evtfile 
A ^ {p', 6' e rotevts : ri < 6' < r2} 

9 -h- max p — min p 

p',S'eA p',5'eA 

if 9>7i: 

9 ^ max ip' + TT mod 27r) — min ip' + vr mod 27r) 
return A, 9 

Table 2: A pseudo-code implementation of the fullAnnuli routine in mathematical shorthand. 

These two pieces of pseudo-code, programmed in Python, provide the ability to make pre- 
cise event regions of multiple radial distances simultaneously. Partial annuli are particularly 
useful for directional analysis, or for selecting very narrow regions for analysis without the 
use of masking files. The full annuli code is particularly useful in central cluster observa- 
tions, where the creation of single complete annulus is possible. An annulus defining system 
is necessary for the use of cluster mass mixing models. 

The rotating algorithm may be described as follows: Let po, Sq be the right ascension and 
declination, respectively, of the center observation. The position of these coordinates on the 
unit sphere is the vector z' = (cospo cos (5o, sinpo cos (5o, sin^o)- As we want this coordinate 
to be rotated to the north pole of the unit sphere, it is desired to rotate the point around 
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a pivot on the equator of the unit sphere. The coordinates of this pivot point are po + 
right ascension, dechnation. The vector coordinates on the unit sphere of the pivot point 
are y' = (— sinpo, cospo, 0). The crossproduct is given by: 



x' = z' X y' 



= (— cos Po sin 60, — sin po sin 60, cos 60) 
Together x', y' and z' form a 3 dimensional rotation matrix R 



x' 



R = \ y' 



z' 



- COS Po sin (5o — sin po sin 5o cos 5o 

— sin Po cos Po 

cos Po cos (5o sin po cos (5o sin 5o 

The action of this matrix rotates po, (5o to the z axis and po + 7r/2, to the y axis. To rotate 
a particular event with coordinates p, (5,and corresponding unit vector v, the dot product R: 

(sin b cos (5o — sin 5o cos 5 cos (p — po) 
cos b sin (p — po) 
cos 5 cos 5o cos (p — po) + sin 5 sin 5o 

It can be seen that p', the rotated right ascension of the event being rotated up to the north 
pole is given by: 



p = arctan 



cos 5 sin (p — Po) 



sin 5 cos 60 — sin 5o cos 6 cos (p — po) 



In the interest of computing annuli sections, it is necessary to compute the great circle 
distance of p, 6 from po, 5o- This is given by the Vincenty great circle distance formula [15] 
which is associated with the new declination of the rotated event: 



S' = arctan 




cos 6q sin (p — po))^ + (cos 6 sin 6q — sin 6 cos 60 cos (p — Po))^ 
sin S sin Sq + cos S cos So cos (p — po) 



Since both p' and S' have many terms in common, they can be computed together efficient. 

The resulting algorithm is implemented in the rot function in the f ilter_annuli .py 
module. It is describe below using pseudo-code in Table |3j 



2.4 Objective and Background Spectra Generation 

Spectra were generated using the ESAS mos-spectra task. CCDs to be excluded were 
automatically loaded as described in §2.2. Instrument and exposure number were also au- 
tomatically generated from processed data. For all observations of Perseus, images were 
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rot{po,So,p,S): 

p' ^ arctan ( cos^sin(p-po) 

6' ^ arctan 



return p', 6' 



sin 5 cos (5o~sin Sq cos (5 cos(p— po) 

y/ (cos 5o sin(p— po))^+(cos 5 sin (5o~sin 5 cos Sq cos(p— po))^ 



sin <5 sin <5o +cos 5 cos 5o cos (p— po ) 



Table 3: A pseudo-code implementation of the rot Python function 

generated along with the spectra for the energy ranges of 400 eV through 1.2 keV and 2.0 
keV through 7.2 keV. Such image generation is optional and time intensive, but for the pur- 
poses of this paper, image generation was necessary to provide redundancy and confirmation 
than correctly formatted annuli sections were generated. The mos-spectra task uses the 
now rotated and redefined event files to produce the annular- shaped spectral regions as well 
as necessary intermediate files for the second task, mos-back, which generates the model 
particle background spectra. 

Ancillary response files (ARF) and response matrix files (RMF) are also created for 
each spectra generated. The ARF files provide the ability to use the spectra generated in 
any standard high-energy spectra analysis software by calculating the effective area of the 
detector in question, while the RMF similarly follows preset standard conventions to describe 
the actual behavior of the detector as a function of energy, through the use of the calibration 
files (plj). 

Object spectra, background spectra, ARF, and RMF files are grouped together as one 
main spectra file through the use of the grppha program. This permits the quick and easy 
loading of all the necessary spectral files into whichever analysis program is to be used, 
simply by specifying the grouped spectra file. 

2.5 Automation Framework 

The tasks listed in the previous sections consume a large portion of time and computer 
power to properly generate. However, the process itself can run into errors in various places. 
The use of the GNU make automation program was indispensable. It permitted this inves- 
tigator to create a single recipe outlining each step that needed to be taken before the next 
task in the pipeline should be run. At the end of each task, a file would be generated should 
no fatal errors occur. The existence of these " completion files" were listed as dependencies on 
the next set of processing tasks. If the file wasn't created, the program would halt analysis at 
the point of failure. This permitted extensive debugging procedures that previously weren't 
available, as well as greatly reducing user error while compiling the parameters necessary 
for each task. This modular design is extremely portable and can be run on nearly any 
UNIX-based computer, greatly simplifying the ESAS procedure. One major benefit of this 
system is the capability to do batch processing for numerous spectra at once, which reduced 
processing time significantly. 
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2.6 Spectra Analysis 



2.6.1 Spectra Preparation 

Spectra were analyzed using tlie 12.7.0 release of the XSPEC spectral analysis software 
package. Local NFWMASS and CLMASS models were used for the purposes of this analysis. 
Cosmological constants in use were Hq = 70, = 0.73, = 0.27, and relative solar 
abundance used for all abundance fitting parameters. Systematic error was fixed at 0. 

Spectra were loaded in groups of three or four as defined by the observation the spectra 
originated from. The spectra analyzed were derived from the full annulus mode of the spectra 



generation tasks (^2.4), but were grouped by observation to preserve directionality. 

Bad event responses were ignored immediately. Any events below the range of .6 keV 
and above 8.0 keV were excluded from fitting, due to the falloff at higher energies and the 
non-standard response of the MOS detectors below ~ 0.4 keV. 



2.6.2 Fitting Procedures 

Fits were run using the standard statistical package loaded into XSPEC; statistics were 
used for the fitting. 

All fits were done after first varying to parameter to be fit until the resultant value 
was reasonably decreased from the default parameter value. In cases where background 
parameters were fit, the program was left to attempt a fit without any adjustment initially, 
due to the large number of parameters needed to be varied simultaneously to achieve a good 
fit (~ 40 parameters). 



2.6.3 Soft Proton, Cosmic Background, and Instrumental effects 

As described in the ESAS Cookbook [14], the background was modeled by multiple 
components for local foreground, instrumentation, and cosmic background sources. Gaussian 
model components were used to describe the two radiation sources (Al K-a and Si K-a), with 
energies of roughly 1.49 and 1.75 keV. For each spectra, a power-law model was added with a 
photon index fixed at 0.2 across all spectra, as recommended. These power-law components 
were not folded into the instrument effective area, depending only upon the RMF file of each 
spectra. 

Cosmic background was modeled through the use of ROSAT All-Sky Survey archival 
data. An annulus of the region between 1 and 1.5 degrees from the center of the cluster was 
added to the loaded spectra and fit directly for the first four channels of the detector [14] . 
This provided the necessary low energy cosmic background filtering. The models applied 
to this spectra included a low energy (< .1 keV) thermal APEC component (unabsorbed), 
an absorbed thermal APEC component of a slightly higher energy (~.25 keV), and an ab- 
sorbed power law fixed at a photon index of 1.46 to account for background cosmological 
X-Ray sources (normalization fixed to the canonical value 8.88 x 10^). Once the background 
model parameters were well fit, the RASS spectrum itself was excluded, as the spectra used 
interfered with the fitting of the mixing models used in the mass analysis. 

Two Gaussian components, centered at 0.65 keV and 0.56 keV were also included to 
account for solar wind-ionized particle contamination, but almost universally these models 
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return extremely low, non-fitting normalizations and were locked at zero intensity. 



2.6.4 APEC and VAPEC Cluster Emission fitting 

Once the background had been well fit, an absorbed thermal APEC component was loaded 
as the primary model for all source spectra. Redshift for these APEC components was fixed 
at the standard valued for Perseus of 0.0183 [T^ [5]. Abundances and temperatures for each 
spectra were allowed to vary freely. The absorbed component was set to values derived from 
the X-Ray Background Tool used to generate the ROSAT background spectrum. 

Once the standard APEC model was reasonably well fit, with a reduced of ~ 1, it was 
replaced with a VAPEC model. A VAPEC model allows more precise control of the abundances 
of 12 different elements. For the purpose of this study, all alpha elements [7j were linked to 
the abundance of oxygen, which was allowed to vary, while Fe was allowed to vary freely. 
All other elemental abundance parameters in the VAPEC model were frozen at zero. Similar 
to the APEC model, all normalizations were frozen to the innermost annuli spectra. 

APEC models were chosen due to the ability to select one central abundance as well as the 
ability to switch between the standard model and the variable abundance model with ease. 
Further, the APEC model is based upon the principle that the source emission, which is an 
ionized plasma contained in the ICM of Perseus, was ionized through collision processes. Such 
a model, based upon thermal bremsstrahlung and line emission processes, is a reasonable 
emission model for the Perseus cluster [2]. 



2.6.5 Mixing Model Fitting 

Following the fitting of the APEC and VAPEC models, the models were further edited to 
include the mixing models NFWMASS [9j and the CLMASS model. 

NFWMASS was initially fit as a simple test to confirm the spectra were properly formatted. 
In one case, the spectra had some slight overlap between annuli sections, resulting in errors. 
This spectra was excluded from all analysis. To provide a good fit for the mass density 
Po/(^(l + ^)^) and total mass of the region Mtotai (Equations ^ and (|3])), temperatures 
were linked to the proper temperature components of the model, while the scale length a 
and potential U were varied manually initially to roughly expected values, to ensure that 
the fitting procedure would converge. The thermal models were absorbed APEC models 



([: 2.6.4), as opposed to VAPEC models, so as to reduce the uncertainty of the large number of 
parameters in the VAPEC model. 

U = inGpoo^fimH (2) 
_ r 47rrVo dr 

^Vltotal — j r Q I rN2 
a \ a' 

Once NFWMASS models were well fit, the NFWMASS component was exchanged for a CLMASS 
component, for the derivation of electron number density and gas densities. Similar tem- 
perature linking procedures were utilized, while densities for each shell were allowed to vary 
freely after being set to roughly expected values [13]. Densities for each region were produced 
in the units of (li]): 

^ ^ keV 
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where u is in units of arcseconds. Conversions to cgs units can be found in ^3.2 



In each case, beta models were not used, as the outer annulus was taken to be the 
approximately outer emission from the cluster. Normalizations for each spectra were tied to 
the normalization of the innermost annuli, as this normalization describes the flux through all 
of the shells of the spherical cluster. The inner radial distance from the center of the cluster 
were in units of arcseconds. To describe the angular extent of each region the AREASCAL 
keyword within the spectra file was edited to the appropriate value for each spectra, while the 
outer edge radial distance of each annuli was stored in the XFLTOOl keyword, as specified 

m. 



3 Results 

3.1 Tabulating Results 

From Table [5} it can be seen that much of the analysis was plagued by high levels of 
uncertainty in the parameter fit, despite generally acceptable statistics. In the case of the 
VAPEC model (Table |4]), this was generally not a significant issue, though it was apparent that 
the linking of so many alpha elements together added an unnecessary level of uncertainty 
into that parameter. Iron abundances tended to be reasonably well fit by comparison, 
with significantly lower uncertainties in almost all regions. CLMASS and NFWMASS model fits 
experience high levels of error (Table [s]). In comparing density, gas fraction, and baryon 
fraction results with Simionescu et al. and Sanders et al., it is apparent that the results 
obtained were generally within a factor of ten of the previous studies' reported results (see 



RA 


DEC 


^1 


R2 


Temperature (keV) 


a 


Error (+/-) 


Fe 


49.927795 


41.234965 


800 


1200 


5.26472 


0.189976 


0.139367 


0.229348 


49.950017 


41.12178 


1200 


1600 


5.29385 


0.121239 


0.100111 


0.170956 


49.948232 


41.014835 


1600 


2000 


4.97118 


4.92575 X 10-5 


0.139413 


0.12039 


49.970342 


40.912062 


2000 


2400 


4.21623 


0.332269 


0.312645 


0.240518 


49.929748 


41.779298 


800 


1200 


6.53401 


0.173715 


0.199214 


0.140577 


49.931043 


41.89284 


1200 


1600 


6.89425 


5.84522E-08 


NO FIT 


0.245601 


50.014582 


42.002226 


1600 


2000 


5.15523 


6.53775E-08 


NO FIT 


0.148863 


49.956595 


42.098751 


2000 


2400 


2.71061 


1.89836 


1.25537 


0.0764963 



Table 4: Sample of VAPEC model temperature, a and Fe abundance fit from ObsID 
0085590201. 



3.2 Calculation 

To generate the final density and mass numbers the annular regions, from the provided 
numbers generated by the NFWMASS and CLMASS models, a few calculations had to be made 
to first. Angular and luminosity distances were calculated to provide a value of roughly .367 
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RA 


DEC 


Ri 


R2 


a 


Error 


U 


Error 


Po 


Error 


49.927795 


41.234965 


800 


1200 


104.222 


444.886 


69.3545 


147.579 


1.58 X 10"*^ 


4.81 X lO"'' 


49.950017 


41.12178 


1200 


1600 


104.222 


444.886 


69.3545 


147.579 


9.79 X 10"^ 


5.21 X lO"-'' 


49.948232 


41.014835 


1600 


2000 


104.222 


444.886 


69.3545 


147.579 


NULL 


NO FIT 


49.970342 


40.912062 


2000 


2400 


104.222 


444.886 


69.3545 


147.579 


1.36 X 10""' 


1.86 X 10"** 


49.929748 


41.779298 


800 


1200 


186.405 


1215.18 


95.7897 


266.187 


2.65 X 10"*^ 


4.32 X 10-** 


49.931043 


41.89284 


1200 


1600 


186.405 


1215.18 


95.7897 


266.187 


9.58 X lO"'^ 


6.19 X 10"** 


50.014582 


42.002226 


1600 


2000 


186.405 


1215.18 


95.7897 


266.187 


9.04 X 10"** 


8.44 X lO-** 


49.956595 


42.098751 


2000 


2400 


186.405 


1215.18 


95.7897 


266.187 


7.47 X 10-s 


8.44 X 10"^ 



Table 5: NFWMASS fit for a and U and CLMASS fit for po from ObsID 0085590201 and ObsID 
0305690301 respectively, a and U are the same for all annuli of each observation, as only 
one value is provided for each spectral group. 



kpc/arcsec and 3 x 10^^ meters This conversion factor was used both in the calculation 
of the gravitating mass density derived from CLMASS as well as gravitational potential from 
NFWMASS. A constant factor of 1.6 x 10^^ in units of ergs/keV was included. The value of 
0.62 is used for the mean mass fraction of a full ionized plasma [T], while 1.6 x 10"^'^ grams 
is used as the mass of a proton. 

The gas density from CLMASS, in cgs units of grams/cm^, is given by the equation: 

1.6 X 10-9 - Po 



Pg 



G(.6- 1.6 X 10-24) . (.367-3 X 10 



2l^2 



The electron number density is given by the equation: 



n, = NJV = ^ ^ X (1 



1.6 X 10-24 ^ tuh-Nh Nh 1.6 X 10-24 
While the total mass density of the region is provided (via NFWMASS): 

1.6 X 10-9 ■ U 

Pt 



G(.6 ■ 1.6 X 10-24) . (a ■ .367 ■ 3 x 102i)2 
Finally, the baryon fraction is provided by: 

_ Pg 



gas 



{pt - Pg) 



In comparison to previous studies of the Perseus cluster, notably the aforementioned 
study by Simionescu and a related study by Sanders et al. [131 [12], it has been determined 
that the error from the anticipated results were roughly of an order of magnitude higher, in 
the case of gas mass density, and an order of magnitude too low for the gas fraction (Fig. 
5|. However, due to the expectations of the sources of error and the extent to which these 
errors could alter these results, these error margins are considered to be acceptable. More 
rigorous practices must be utilized in future study of this region. 
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Figure 5: The dependence of the deprojected fgas, the baryon gas fraction, as a function of 
radius. High levels of error skew this results an order of magnitude too low. 



In the case of temperature profiling and iron abundance, via the VAPEC model, results were 
far closer to presently accepted values to within a respectable margin of error [2]. However, 
it is desirable to produce higher quality results for across the entire region, as the size of 
the annular sections used in this study are of significant breadth compared to the whole of 
the observed clustered. This provides a major limitation on the ability to pinpoint exact 
regions and define temperatures, as the temperatures and elemental abundances produced 
by the VAPEC model are not provided as a direct function of radius. Increasing the number of 
annular regions, as mentioned in Q is expected to provide a continuous function describing 
the temperature and abundances dependent on direction and radius as well as significantly 
improved statistics. One significant area of concern, as mentioned below in 3.3[ is the 
distinctly non-uniform temperature distribution of the cluster (Fig. [6]). This, tied with the 
spread of abundances observed (Fig. [?]), provides reason to argue that the Perseus cluster 
does not model a system in hydrostatic equilibrium [6l|5]. However, it is expect that further 
study (Q should yield results with more reasonable error to combat this potential failing of 
the mass mixing models. 

As described in [12], non-spherically symmetric bubbles and ripples in the immediate 
central region of the cluster exist, while in some cases shock waves extending out to ~170 
kpc. This could particularly be indicative of recent interactions and collisions with galaxies 
or the infall of large quantities of matter into the central AGN. Furthermore, there is recent 
evidence of significant gas motion within the immediate 100 kpc of the central AGN [1]. 
With these considerations, it is possible that some error in the model fitting stems from the 
failure of the initial assumption concerning the equation of state of the cluster itself. 
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Figure 6: A profile of the projected temperatures across the Perseus cluster. Circle size is 
indicative of the level of error of the parameter value. Cluster center is at the origin. The 
code referenced in §2.3.1 was used to provide the coordinate rotation to produce both Fig. 
|6]as well as Fig. [7| North is in the direction of the upper left quadrant. 



3.3 Error Considerations 

While the error observed in this study was seen to be unusually high, much of the 
error, it is believed, can be reduced in future analysis. Only one detector was used in this 
analysis, though spectra were produced for both the MOSl and M0S2 detectors initially. 
This current work on the MOSl detector was done simply as a test case to ensure that 
results obtained from this new method were not significantly different, to the point where 
this method of analysis should be abandoned. The elimination of the second detector poses 
unique problems. Statistics were greatly reduced by only using half of the data, averaging 
results across the two detectors was not an option, and the inadequate performance of the 
primary detector due to physical damage. 

The detector used in this study was damaged in a micrometeorite impact, resulting in 
the total loss of the CCD6 of the MOSl detector. In normal operation, such a loss would 
not be of significant detrimental effect to the quality of observation data; annular sections 
could be lacking significant amounts of the effect light collecting area (Fig|9]). This is obvious 
when the statistical robustness of the Fe abundance and temperature fitting from the VAPEC 
models is considered, as this model doesn't rely upon the annular shape for fitting. Large 
gaps in the annular sections are expected to have been a major source of error for the 
CLMASS and NFWMASS models, as mixing models are dependent upon having a continuous set 
of annuli from a defined inner radius to the edge of the cluster (or to the farthest radial 
extend possible). 
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Figure 7: A profile of the Fe abundance across the Perseus cluster. Circle size is indicative 
of the level of error of the parameter value. 
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Figure 8: A profile of the electron number density across the Perseus cluster. In combination 
with Fig. [6] and Fig. [7] are the main examples of high error stemming from the NFWMASS and 
CLMASS model fits. 

Another potential source of induced error was the chosen spacing and angular size of the 
annular regions. For the purpose of producing spectra that would be guaranteed to have 
adequate numbers of events, large regions, some annuli sections with > 30,000 events were 
used. The parameters for the annular sections were set to have radial extent of 400", while 
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Figure 9: On left: Image from the annuli across 1200" through 1600" with damaged CCD 
region in upper right corner. On right: Image from annuli across 1200" through 1600" for 
the intact M0S2 detector. Exposure retrieved from ObsID 0305690301. 



the angular size varied between | to roughly 27r radians. In the course of the analysis, it was 
determined that having too few spectra for each directionality strongly limited the ability of 
the mass mixing models to general accurate fits. Simply, there werent enough data points 
to produce reliable statistics for the gravitational potential and mass density parameters. 



4 Future Study 



The pn-CCD has yet to be tested in this method, though in previous study, its ability to 
adequately deal with background noise was much more limited than the MOS instruments. 
However, in the effort to improve statistical quality of results, its use in this form of analysis 
will be revisited. 



Using the partial annulus mode described in ^ 2.3.1 the angular extent of the spectrum 
derived can be shrunk considerably, while the radial size can also be reduced to at least 
a factor of 4 smaller. Considering just the reduction in radial size, this would potentially 
produce NFWMASS and CLMASS fits with significantly reduced error, providing potentially 4 
or more times the amount of data points for the mixing models to fit. It is also possible to 
extend farther into the core of the cluster, where the quality of the X-Ray data is significantly 
higher, due to both the higher incident flux from this region, as well as the existence of two 
separate observations of the same region. This will permit the ability to greatly improve the 
statistics of the model fittings for the innermost regions of the cluster. 

While this study has been unable to produce results of great statistical value, it has been 
shown that, with some adjustments, this method of analysis could provide great insight into 
the composition, temperature, and matter distributions of "local" clusters. With this in 
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mind, it is logical that the next step is to examine other clusters with a similar methodology. 

Furthermore, the automated procedure outlined in this paper was designed for the purpose 
of wide scale application. Study of clusters such as Coma, Virgo, and Ophiucus are expected 
to take place in the not-too-distant future. 
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